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ABSTRACT 

This thesis investigated a laboratory synthetic aperture sonar designed to test the 
algorithms and techniques needed to detect, classify and identify minelike objects. Previous 
synthetic aperture sonar work at NPS achieved 5 cm range resolution and 1 cm azimuth 
resolution. This thesis developed a pulsed, frequency modulated, synthetic aperture sonar that 
achieved range and azimuth resolutions of about 1 cm. The processed images clearly reveal 
targets with a high degree of certainty. However, the ability to classify and identify mines and 
rocks is less certain because of speckle and glint effects. The high resolution algorithms 
improved the detection and overall image quality of targets, and achieved a signal to noise 
ratio of 35 dB. The 2: 1 frequency spread of the FM chirp increased the signal to noise ratio 
by 20 dB compared to an unfocused synthetic aperture system. However, a significant finding 
is that resolution alone is not sufficient to classify and identify minelike targets in complex 
backgrounds. Resolution of this problem will require a different approach such as utilizing 
adaptive acoustic daylight to avoid the speckle and glint problems inherent with coherent 
illumination. To achieve a classification and identification capability, a completely different 
approach to acoustic illumination and signal processing is needed. 
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I. INTRODUCTION 



The use of synthetic aperture techniques in radar is well established. Such techniques 
have been used extensively in both military and civilian applications starting around 1955. 
Since 1961, extensive literature on synthetic aperture techniques has been published. 
However, the application of synthetic aperture techniques in sonar has been rare. 

There are a number of applications where high resolution underwater imaging is of 
importance. One such application is in minehunting operations where the ability to detect, 
classify and identify a mine using a sonar system provides for safer and faster removal of the 
mine threat. This is the primary motivation for this study. 

This thesis looks at the application of synthetic aperture techniques to achieve a high- 
resolution imaging sonar. It explores the improvement in the quality of sonar images by the 
application of relatively simple synthetic aperture techniques. In addition to improving the 
resolution of synthetic aperture devices, improvement of signal-to-noise ratio is also 
important, and is studied in this thesis. 

The objectives of the study are to analyze and compute the potential improvements 
in imaging quality using selected, synthetic-aperture data processing techniques and to 
understand any physical limitations that hinder wider applications of such techniques in 
sonar. 

The study was carried out with an experimental data collection system that used 
sound propagation in air in the 16 kHz to 32 kHz frequency range. These data were then 
processed using unfocused and focused synthetic aperture algorithms. Previously, an azimuth 
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resolution of 1 cm was achieved while the range resolution was about 5 cm [Welter, 1995]. 
A linear frequency modulation technique was applied in this thesis to attain a range 
resolution of 1 cm, which when applied together with the 1 cm azimuth resolution focused 
algorithm improved the overall image quality. 

The experimental data collection and processing provided a ‘hands-on’ 
understanding of applying synthetic aperture techniques to achieve the theoretical imaging 
resolution. This provided valuable insights into the intricacies of deploying synthetic aperture 
sonar systems. 
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II. BACKGROUND 



A. MOTIVATION FOR THE STUDY 

There are a number of applications where fine resolution underwater imaging is 
critical. One such application is in minehunting. The ability to detect, classify and identify 
a mine using a sonar system augments the safe and fast removal of the mine threat. However, 
there are physical limitations on wavelength and transducer size that suggest synthetic 
aperture techniques should improve the image quality. One such limitation is that a single 
sonar transducer cannot be as big as needed since it is deployed on small ships such as 
minehunters or remotely-operated vehicles controlled by minehunters. A potential answer 
to this problem is a synthetic aperture system since it can provide very fine along-track 
resolution, independent of range and frequency, with a small transducer. 

B. CHAPTER OVERVIEW 

The primary objective of any sonar system is to be able to detect, classify and identify 
contacts. To achieve this objective, a minehunter operator needs high resolution sonar images 
that discern between mines and non-mine objects. Two general categories of processing 
techniques exist. One uses real aperture processing while the other uses synthetic aperture 
processing. 

The key characteristics in comparing synthetic aperture sonar systems and 
conventional real aperture systems such as the side-looking sonar are the image quality and 
the signal-to-noise ratio. Good measures of image quality are the range and azimuth 
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resolutions. This chapter provides an overview of the basic operating principles of real 
aperture sonar, describing the beam shape and pattern, and the range and azimuth resolutions. 
Pulse compression provides better range resolution while synthetic aperture techniques 
improve upon the azimuth resolution. In addition, both the pulse compression and synthetic 
aperture techniques enhance the signal-to-noise ratio. The chapter concludes by looking at 
two interesting phenomena, speckle and glint, which have significant impact on the image 
quality. The effect of synthetic aperture processing and pulse compression on these 
phenomena is studied in this thesis. 



C. REAL-APERTURE SONAR 

Conventional sidescan sonar systems use real-aperture image processing that provides 
two-dimensional reflectance maps of the acoustic backscatter energy from the objects in the 
environment. The system typically uses an acoustic projector that moves in a straight line. 
As the projector transmits acoustic pulses, a hydrophone listens for the echos. The time delay 
of each echo provides a measure of the range while the ping-to-ping motion of the projector 
gives the along-track or azimuth image dimension. 
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1. Range Processing 



The most common operating mode for sonar or radar systems is to transmit 
pulses which will scatter off targets and be detected by a receiver. The pulses are usually of 
very short duration, microseconds or milliseconds. This duration is known as the pulse 
width. The pulses are separated by a longer period, known as the pulse repetition interval 
(PRI), which depends on the maximum desired detection range, as will be explained 
subsequently. 

The strength of the signal from a target is measured by the power scattered 
off the target and received by the system. The received power is 



P = 

r 



P ( G ‘oA 
(47 Ur 2 ) 2 



( 1 ) 



where P r = power collected by the receiver with cross sectional area A, 

P t = power that is sent out by the transmitter with cross sectional area A, 

G‘ = antenna gain, 

ct = backscatter cross-section of the target and 
r = range of the target. 

The denominator factor in brackets represents the spherical divergence of the transmitted 
waves followed by the divergence of the scattered waves. 

Pulsed systems determine the range by measuring the transmit/receive delay 
interval. Figure (1) depicts both the time line for the transmission of a pulse of pulse width 
x with a pulse repetition interval T and the echos from two closely spaced objects. 
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Figure 1 Time Domain Representation of a Transmitted Pulse and 
Corresponding Echoes 



The two objects separated by range 5R are resolved if their echos do not 
overlap in time as shown in Figure (1). The round trip time for an object at a range R is given 
by 

2R 

t = — , ( 2 ) 

c 

where c is the propagation speed in the medium. 

A second object that is 8R in range from the first object introduces an 
incremental delay 5t given by 



, * 8 , . 2 (*+ 5 *> 

c 



( 3 ) 
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Subtracting Equation (2) from Equation (3) gives the temporal separation of 



the two echos as 

{, = M. 

c 



(4) 



A point target will generate an echo that is x in duration. Echos from two 
point targets must be separated in time by at least x to be distinguishable from one another. 
Therefore, we set 5t equal to t. The range resolution of the pulsed system is the minimum 
distance that two objects must be separated in order to be discerned. This is obtained simply 
by solving Equation (4) for 8R 



6 *« 



ct 

2 



(5) 



Equation (5) indicates that the range resolution is proportional to the pulse 
duration, which means that shorter pulses will provide finer resolution. However, the 
transmission of shorter pulses means that less energy is emitted, which reduces the detection 
ranges. Therefore, in a real aperture sonar, it is not possible to have very fine range resolution 
at long ranges without employing pulse compression to reduce the effective length of the 
pulse. 

Figure (1) also shows that the PRI has to be longer than the desired maximum 
range in order to avoid ambiguity. The PRI has to be set sufficiently large so that all the 
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echos from a pulse are received prior to transmitting the next pulse to circumvent the 
problem of range ambiguity. This is equivalent to requiring that the time between pulse 
transmissions be greater than the time for the returns from the far-edge, or maximum 
processed range, of the beam footprint. 

The minimum acceptable PRI, and hence the maximum pulse repetition 
frequency (PRF), to avoid ambiguous range returns are 



PRI in = 

mm 



2 R 



( 6 ) 



and 



PRF 



max 



PRI . 

mm 



( 7 ) 



respectively. 

2. Pulse Compression 

One way to overcome the limitation of having range resolution dependent on 
short pulses with low energy as suggested by Equation (5) is to adopt a pulse compression 
technique. This entails sending out long coded pulses that are either frequency or phase 
modulated. This technique is also commonly known as chirping the pulse. The compression 
process correlates the echo pulse with a delayed copy of the transmitted waveform as shown 
in Figure (2). 
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Figure 2 Pulse Compression Principle 

This is equivalent to separating the echo waveform into parts based on the modulation and 
delaying each separate part such that they are made simultaneous in time. The component 
parts, which occur in the time sequence in the echo, are summed at the same time producing 
a narrower pulse of higher amplitude. For a signal with m data elements and a transmitted 
waveform of n elements, this correlation will result in a compressed signal with (m+n-1) data 
elements. 

The linear frequency modulated waveform consists of a rectangular transmit 
pulse of duration T as shown in Figure (3a). 
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Figure 3 Linear FM Waveform and Processing 



The carrier frequency is swept linearly or chirped over the pulse length by an amount Af. The 



frequency, f(t), of the carrier f 0 will be modulated according to 



AO = 4 - ^ 




( 8 ) 



As shown in Figure (3b), the higher frequency components of the received 



signal experience longer delay times than the lower frequency components. The received 



signal s(t) will have the form 
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( 9 ) 



5(r)=cos[2wft--^)f] M<- . 

2 T 2 



where the argument of the cosine function is the phase of the received signal. 

Eaves & Reedy [Eaves & Reedy, 1987] show that the matched-filter impulse 
response for the signal s(t) shown in Equation (9) is a time inversion of s(t) 



which is a sine function with a phase factor. 

The output signal vj/(t)ffom the pulse compression may be characterized by 
an envelop having a higher amplitude and a narrower pulse length than the transmitted 
envelope as shown in Figure (3c). Eaves & Reedy show that the compressed pulse length, 
x c , and the compression ratio (CR) are 




( 10 ) 



Furthermore, Eaves & Reedy show that the correlation of s(t) and h(t) is 




( 11 ) 



T 




( 12 ) 



c 
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and 



• ( 13 ) 

C 

In summary, pulse compression enables use of the energy of the long 
transmitted pulse for the purpose of detection and, at the same time, the bandwidth of the 
short compressed pulse for range resolution. The range resolution attained using pulse 
compression is 



6 R 



C 




C 

2 BW 

c 



(14) 



where x c is the compressed pulse width and BW Represents the bandwidth of the compressed 
pulse. 

Besides better range resolution, pulse compression processing improves the 
signal to noise ratio. On a single-pulse basis, the probability of detection at a certain range 
increases with the signal-to-noise ratio. A pulse compression system achieves higher signal 
to noise ratios from the increased energy at no expense to range resolution by utilizing a 
longer transmit pulse. 

3. Azimuth Processing 

For processing in the azimuth or along track axis, it is necessary to look at the 
beam pattern of the transducer for both transmission and reception. The beam pattern is 
determined by the transducer’s size, shape and frequency. One common measurement of the 
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beam pattern is the 3 dB beamwidth. This is the angular width of the main beam where the 
response is within the 3 dB or half-power region of the peak response. 

The half-power angular beamwidth 0 3R of a uniformly-weighted, rectangular 
aperture of length D or circular aperture of diameter D depends on the 2 dimensional Fourier 
Transform of the field across the aperture, which gives 

Q 3 R =X/D ( radians ) , (15) 

where X is the acoustic wavelength [Edde, 1993], 

Given the 3-dB beamwidth, a target within this region with sufficient target 
strength will be detected. From the Pythagorean theorem, the distance to a single stationary 
target as measured by the sonar as the aperture moves along a straight track, will be a 
hyperbola as shown in Figure (4). 
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Azimuth 

Figure 4 Detection of a Single Stationary Target 
The target will be detected at varying ranges as the aperture moves pass it. The detection 
range, R, of the target at the indicated azimuth position is 

R = + A? , (16) 

where = range of target at closest point of approach and 

Az = azimuthal distance from the point of closest approach. 

The target data is digitized at uniform time steps and stored as a matrix where 
each azimuth position represents a column and the range bins are rows. The azimuth/range 
matrix is schematically shown as grids in Figure (4). For the single stationary target shown 
in Figure (4), the transducer will see returns for about seven azimuth positions. Figure (4) 
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only displays detection of this target at three of the azimuth positions and the corresponding 
detection ranges. 

In real-aperture sonar, the actual beamwidth at range R is a measure of the 
along-track or azimuth resolution. The criterion for azimuth resolution is that objects at the 
same range separated azimuthally by more than the beamwidth are resolved [Edde, 1993]. 
Thus, the azimuth resolution for a real-aperture sonar is 






RX 

D 



(17) 



Equation (17) shows that real-aperture systems require very narrow beam widths in order to 
achieve fine along-track resolution. To keep 5A small as range increases, the frequency of 
the system and/or the physical aperture length must increase. However, even with a narrow 
beamwidth, the azimuth resolution is still proportional to the range. This is where synthetic 
aperture sonar overcomes both the physical limitation on how big the aperture can be and 
provides an azimuth resolution that does not vary with range. 



15 



D. SYNTHETIC APERTURE SONAR (SAS) 



The primary motivation for using SAS centers on its ability to generate very fine 
azimuth resolution that is independent of both range and frequency, without the use of a large 
transducer. In fact, for SAS, the smaller the size of the transducer, the better the azimuth 
resolution, which is just the opposite behavior found for a conventional sonar. Range 
processing and pulse compression are the same as that for real aperture sonar as explained 
above. 

SAS is based on the generation of an effective long antenna by adding the signal 
amplitude coherently rather than by the use of a big transducer. In other words, synthetic 
aperture processing relies on increased signal processing complexity as opposed to large 
physical apertures to obtain fine azimuth resolution. 

The physical aperture of a SAS may be regarded as one element of a linear array 
extending in the direction of platform motion as shown in Figure (5). 
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where 0 3S is the 3 dB beam width of the SAS. The 3 dB beam width for a synthetic aperture 



is 

A A 

8 » = Tl ’ <»*> 

where L is the length of the coherent addition and the factor of two arises from measuring 
the round trip phase of the returned signal. 

The azimuth resolution of the synthetic antenna at range R is then 



6 A 



s 



xr 

2 L 



(19) 



The signal processing approaches for SAS are categorized generally as unfocused and 
focused apertures depending on the treatment of the phase variation of the signals received 
by the transducer. 

1. Unfocused SAS 

Figure (6) shows the formation of a synthetic aperture with the sensor moving 
along a nominally straight line with constant velocity for a distance L ; , referred to as the 
coherent integration distance. Averaging of consecutive sonar returns without compensating 
for the spherical waveform will provide an azimuth resolution intermediate between that of 
conventional real aperture beamforming and that of focused synthetic aperture beamforming 
[Heering,1984], This situation is equivalent to assuming that the target is in the far field of 
the aperture. 
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Azimuth 




Range 

Figure 6 Synthetic Aperture Formation Geometry 

The phase associated with each echo changes with time since it is proportional 
to the round-trip distance. Thus, a phase variation proportional to 8R„ for the n* element 
exists either side of the point of closest approach. An unfocused aperture results by ignoring 
the phase variation in the data collected over the aperture path. One assumption made for 
unfocused synthetic aperture system is that the maximum round-trip phase variation should 
not exceed n/2 or X/8 for the duration of the quasi-coherent summation [Bruce, 1992, Eaves 
& Reedy, 1987, Curlander & McDonough, 1991]. Another assumption states that this phase 
variation should not exceed n/4 or X/\6 instead [Edde, 1 993] . In both approaches, assuming 
is small, the maximum length of the synthetic antenna for coherent integration will be 



For 8 R < — , L 2 < XR , 
" 8 



(20a) 
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and 

p X I XR 

For bR * < 77’ 1 * “T - ' ( 20b > 

1 0 Z 

With the length of the synthetic antenna limited by either Equations (20a) or 
(20b), utilizing Equation (18), the azimuth resolutions for unfocused synthetic aperture 
system are 



For 6 R 



5 A, 



V^X 

2 



(21a) 



and 



For 6 R < — , S.4 = A 

" 16 s \ 



Rl 

2 



(21b) 



As stated above, formation of an unfocused synthetic aperture image involves 
the coherent summing of returns collected over its synthetic length L. Edde also suggests that 
a moving sum of a chosen number of azimuth elements at a fixed range is also a form of 
unfocused SAS imaging. The process involved in the moving sum is illustrated in Figure (7). 
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Figure 7 Unfocused SAS Processing 

The algorithm stores echo data for each transducer position up to the maximum range of 
interest, which creates a column of data. The number of columns or transducer positions to 
be summed will be based on the coherent integration length or maximum phase variation 
adopted. After summation over the coherent integration length, each summation is stored as 
a column in the display array. This display matrix gives one unfocused synthetic aperture 
image. 

2. Focused SAS 

When the coherent aperture length is sufficiently long so that Equation (20a) 
or Equation (20b) no longer holds, the phase variations from the hyperbolic range variation 
must be removed. This is equivalent to the Fresnel or near field situation in optics. Focusing 
of the synthetic array corrects the phases of the received signals so that all the elements of 
the resulting array are at the same distance or phase from the source. This is illustrated in 
Figure (8) which shows the returns from a point target. 
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s 

Azimuth 

Figure 8 Focusing for a Point Target 

Focusing of the point target at a particular combination of azimuth and range 
cells involves the correction of the phases and summation of the shaded elements as shown 
in Figure (8). Applying the Pythagorean theorem, the phase correction for the n* array 
element for the image cell shown in Figure (8) is 

(6 R n + R) 2 = R 2 + ( nS f , (22) 

where R = range from the broadside SAS element to the position being corrected, 

8R,, = range difference between broadside SAS element and the nth element, 
n = the number of elements being corrected and 
S = spacing between elements. 
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Solving Equation (22) for 5R„, assuming (5R„/2R) « 1 , gives 






n 2 S 2 
2 R 



( 23 ) 



where the two in Equation (23) accounts for the round trip for the wave propagation. 

Figure (9) shows one step in recovering the image [Edde, 1992]. The data 
array shows the signal collected for each range bin (rows) and for each azimuth bin 
(columns).The mask represents the real antenna beam and the data within the mask is the 
data for one S AS coherent summation. Phase corrections are applied to the entire data array 
within the mask using the method explained in Figure (8) . The resultant array summation 
extends across the data points along the hyperbolic curve. This sum produces one image 
pixel at the range and azimuth being processed. The data array is then stepped one image 
pixel distance further along the azimuth path of the sonar and the above process repeated for 
the new combination of range and azimuth cells. After producing an entire row of azimuth 
pixels, the process must be repeated for the next range bin and the entire set of azimuth bins. 
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Figure 9 Focusing Process for One Azimuth Element 
Effectively, a focused array coherently adds the signal amplitudes from all the data elements 
along a curved hyperbolic arc. 

The maximum synthetic aperture length is limited by the 3 dB azimuth 
beamwidth at a given range of the actual transducer. This is because targets that are outside 
the main beam lobe will have relatively smaller returns. The maximum synthetic aperture 
length is therefore determined by the footprint on the ground and is 



L 



max 



RQ 3S 



xr 

D 



(24) 



Substituting L max into Equation (19) gives the resolution of the sonar as 



6 A s = DU . 



(25) 
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Equation (25) is a remarkable result that explains the interest in implementing 
a focused SAS algorithm. It shows that the focused SAS azimuth resolution is independent 
of both wavelength and range. As range increases, so does the maximum length of the 
synthetic aperture. It is also independent of wavelength because the actual beamwidth 
illuminated by the transducer A/D increases with longer wavelength as does the effective 
aperture length L. In addition, the azimuth resolution gets better with a smaller real aperture, 
exactly the opposite behavior seen for a real antenna. The A/D broadening of the beamwidth 
with a smaller physical aperture is counteracted by the larger illuminated region on the 
ground, which increases the length of the synthetic aperture. There is, of course, a limit as 
to how small the real transducer can be since it must provide enough gain and aperture to 
assure an adequate signal-to-noise ratio. 

3. Signal-to-Noise Improvement Factor 

While range and azimuth resolutions are important determinants for image 
quality, an equally important aspect is the signal-to-noise ratio. In sonar, the signal-to-noise 
ratio is usually given in decibels. A synthetic aperture sonar coherently adds a number of 
pulses, namely, the pulses occurring during the time necessary for the sonar to move a 
distance equal to the synthetic aperture length. The signal-to-noise improvement factor due 
to synthetic aperture generation is proportional to the number of pulses integrated [Cutrona]. 
In addition, the improvement factor from pulse compression is proportional to the ratio of 
uncompressed pulse length to the compressed pulse length. For a system adopting both 
synthetic aperture and pulse compression processing, Cutrona expresses the overall 
improvement factor in dB as 
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IF = lOlog 



( 26 ) 






10 

V 




PRF_RX] 

2b A V J 



where x, and x 0 are the uncompressed pulse length (in seconds) and the compressed pulse 
length respectively, V is the speed of translation of the sonar aperture (in feet/second), R is 
the range (in feet) and X is the wavelength (in feet). The second term in brackets in Equation 
(26) is the number of pulses integrated. 



E. SPECKLE AND GLINT 

In images obtained from sonar and radar systems, the effects of speckle and glint are 
significant. These two phenomena are also of particular interest to this thesis because they 
become predominant as the target size gets small and the resolution improves. 

L Speckle 

When nearly monochromatic energy is reflected from an object with a rough 
surface, the wave resulting at any moderately distant point consists of many coherent 
components or wavelets, each arising from a different microscopic element of the rough 
surface. 

Interference of the dephased but coherent wavelets results in the granular 
pattern of intensity that is called speckle as shown in Figure (10). The range variations create 
localized destructive and constructive interference which appears in the image as bright and 
dark speckles or granulation. As a result, the detailed structure of this granularity bears no 
obvious relationship to the macroscopic property of the illuminated object. 
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Figure 10 Images of a rough object formed with coherent light 

2. Glint 

Many targets of practical interest are complex targets consisting of many 
scatterers. The backscatter energy from a complex target is a function of how the 
contributions from the scatterers on the target recombine at the transducer. Glint is the result 
of small local target components that act like comer cubes that reflect a relatively larger 
amount of the incident radiation back to the source. 

In the case of a single point target, like a sphere, the apparent location of the 
target is the center of the spherical wavefronts re-radiated by the target. If the target is 
complex, the wavefronts are no longer centered on a single point. 
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Besides the effect on target position, glint will also obscure details of the 
target in the neighborhood of the target points with strong returns. The effect of glint is 
especially prevalent in the case of a complex target that has scatterer points of widely 
different strengths. The strong return and its associated sidelobes can obscure subtle details 
of the target. 
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III. SYNTHETIC APERTURE SONAR EXPERIMENT 



A. HARDWARE 

Figure (11) shows the schematic setup for the experiment [Welter, 1995]. 




Figure 11 Equipment Setup for Experiment 

The synthetic aperture sonar was developed by Prof D.L. Walters and described in 
previous thesis work [Welter, 1995]. It consisted of a linear track 1.2 meters long that 
includes a small trolley pulled by a cable. The trolley carried a 1.3 cm electrodynamic 
transducer from a Sony earphone and a 0.25 cm Larson Davis microphone. During operation, 
a DC motor pulled the trolley at a rate of about 1 cm/sec from one end of the track to the 
other. 

The transducer was driven by a HP 33 120A function generator that was configured 
to produce a sequence of sine wave bursts or a frequency modulated chirp. The microphone 



29 



included a Larson Davis 2200C preamplifier followed by a Stanford Research SR560 
amplifier and a Wavetek 735 band pass filter. The filtered output drove an AD AC 5508 high 
speed, 1 2 bit, A-D converter contained within an IBM PC, 486-66 computer system. 

During use, crosstalk between the transducer and the microphone was sufficiently 
high that it was necessary to introduce a 2 cm by 10 cm wood baffle between the source and 
receiver to suppress the direct interaction. 

The data acquisition process involved moving the transducer and receiver on the 
trolley. The objects were placed within one meter from the track and about mid-way on the 
track. As the trolley moved down the track, the object was illuminated at a PRF of 1 8-20 Hz. 
The sampling rate was 1 00 kHz and each sample was 12 bits long. A total of 584 x 1024 data 
elements were captured for each run for a total of 1 .2 MB per image. The time for each run 
was determined using a stop watch and was typically 60 to 70 seconds. The along track 
distance traveled by the trolley was measured using a ruler mounted on the track. 

B. SOFTWARE 

The data collection software was a commercial product, Testpoint version 1.1. 
Appendix A contains the program used for the data collection. This program consists of a 
loop of 1024 azimuth positions, each of which collects an array of 584, 12 bit data elements 
at a rate of 1 00 kHz. Each 584 data element for a transmitted pulse was collected and stored 
to disc before the transmission of the next pulse. 

The synthetic aperture algorithms were written in Matlab version 4.2c. 1 , which are 
contained in Appendix B. 
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C. KEY PARAMETERS & PERFORMANCE 

Initially, and in previous work [Welter, 1995], the function generator was 
programmed to send out bursts of sinusoidal waves to the transducer. This was changed by 
downloading a FM chirp to the HP 33 120A function generator as shown in Appendix C. The 
transducer had a nearly uniform frequency response in the 1 5 kHz to 40 kHz frequency 
range. A uniform response would provide the best pulse compression of the collected data. 
The other key parameters selected for the experiment are explained below. 

1. Fixed Frequency 

The fixed frequency pulse transmissions were at 25 kHz, which had a 
wavelength of 1.37 cm. A burst of five cycles provided good range resolution of 3.4 cm and 
good signal-to-noise returns. The pulsewidth of each transmission burst was therefore 0.2 ms. 
For a one meter maximum range and a speed of sound of 343 m/s, the interval between 
pulses could not be less than 5.83 ms. The burst rate was set at 1 8 per second which meant 
that the time between bursts was about 55 ms. 

The diameter of the transducer was about 1 .3 cm. With a frequency of 25 kHz, 
the 3 dB beamwidth was about one radian. This was found to be consistent with 
measurements carried out in an anechoic chamber [Welter, 1995]. 

The Nyquist sampling criterion states that in order for information to be fully 
recoverable after sampling, it must be sampled more than twice for each sinusoidal cycle of 
information present. The sampling rate must be greater than twice the highest frequency 
present. In this experiment, the sampling rate of 100 kHz was four times that of the 
transmitted frequency, which more than satisfied the Nyquist criterion. 
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Given the above parameters, the data array collected for one run was for 584 
range bins and 1024 azimuth bins. The range resolution was range independent. For a five 
cycle sinusoidal burst, the pulsewidth was 0.2 ms and range resolution was 3.4 cm. 

For the unfocused SAS processing, using Equation (20a), the azimuthal 
integration length could not exceed 12 cm at 1 m range and 8 cm at a range of 0.5 m. The 
azimuth resolution for real-aperture processing was range dependent and was 105 cm at lm 
and 53 cm at 0.5 m. The unfocused SAS algorithm gave a range dependent azimuth 
resolution of 5.9 cm at 1 m and 4. 1 cm at 0.5 m. The focused SAS algorithm gave an azimuth 
resolution of 0.7 cm that was range independent. 

2. Chirp 

The pulse compression waveform used an up-chirp frequency band that 
ranged from 16 kHz to 32 kHz that had a compressed pulse bandwidth of 16 kHz. The 3 dB 
beamwidth was narrower in this case, about 22°. The data sampling rate was at 100 kHz 
which still met the Nyquist criterion. The data array for one run was also 584 range bins and 
1024 azimuth bins. The azimuth resolutions were the same as for the fixed frequency case 
since pulse compression only alters the range resolution. With pulse compression, the range 
resolution was 1.1 cm (0.8X) compared to 3.4 cm for the uncompressed pulse. 
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IV. DATA PROCESSING & RESULTS 



For the tests, several objects were used to evaluate the processing algorithms. These 
objects included a sphere, batteries arranged in a shape of an N, rocks, an anti-invasion mine 
and anti-personnel mine. The sphere provided a point target for evaluating the range and 
azimuth resolutions, and the signal-to-noise ratios. This was the primary test target used to 
fine-tune the processing algorithms since the expected image was a point. 

After developing the data processing algorithm on the sphere target, it was possible 
to proceed to more complex targets. The N-target provided a variable range and azimuth 
geometry needed to evaluate the range and azimuth resolutions. The rocks, anti-invasion 
mine and anti-personnel mine represented severe challenges to the data processing 
algorithms in its ability to detect, classify and identify these targets. This last group of targets 
could be encountered operationally and were of special interest in this study. 

A. POINT TARGET - SPHERE 

The sphere, 12 cm in diameter, was a simple point target that was useful to 
evaluate the performance of the processing algorithms in terms of signal-to-noise 
enhancements, as well as range and azimuth resolutions. A picture of the sphere with a 15 
cm ruler in the foreground is shown in Figure (12). 
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Figure 12 Sphere with a 15 cm ruler in the foreground 
Such a known target augments the development of the processing algorithms 
with steps that are not obvious from the literature on SAS processing. The steps that were 
found necessary were the correction for a DC offset in the A-D converter, range amplitude 
compensation, selection of the FM chirped signal waveform and the choice of the translation 
distance. These are explained in detail below. 

1. Data Offset 

The data collected using the ADAC 5508 analog to digital converter was 
found to have a DC offset of about 4 mV and this was enhanced by the range dependent 
signal amplification process as discussed below. The processed signal, with signal 
amplification, in the range axis is shown in Figure (13a). This clearly shows a range 
dependent bias. Figure (13b) shows the same data after the offset was subtracted. For all 
subsequent processing, the DC offset correction was applied. 
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13a Focused SAS Data without a DC Offset Correction 
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13b Focused SAS Data with a DC Offset Correction 

Figure 13 Focused SAS Data of Sphere with and without a DC Offset 

Correction 
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2. Signal Amplification 



The acoustic pressure declined by r' 1 (1 /range) from spherical spreading 
traveling from the source to the target and had another r' 1 degradation on the return path. The 
total round trip loss was proportional to r' 2 . The coherent, synthetic aperture processing of 
the received signal alters this behavior. In essence, the receiver is a phase array where N 
number of pulses add coherently. This number of elements included in the coherent 
summing is proportional to range because of the XfD spreading of the energy. This 
compensates for the return path r' 1 spreading of the acoustic pressure wave and the total range 
dependence of a focused, synthetic aperture sonar image is r' 1 . The unfocused synthetic 
aperture range dependence is intermediate between r' 1 and r‘ 2 . 

Figures (14a) and (14b) show the unfocused SAS processed signal in the 
range axis without range correction and the same data with range correction (r 2 ) respectively. 
The high return at the near range is from the cross talk between the transducer and the 
microphone. The effect of the range correction reduces this direct signal input. The 
unfocused SAS images of the sphere shown in Figure (15a) and (15b) illustrate this more 
clearly. In the following figures, the 584 range bins have been block-averaged in groups of 
8 to produce 73 range bins and the 1024 azimuth bins have been block-averaged in groups 
of 8 to produce 128 azimuth bins. 
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14a Unfocused SAS Data without Range Correction (r 2 ) 
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14b Unfocused SAS Data with Range Correction (r 2 ) 

Figure 14 Unfocused SAS Data of Sphere with and without Range 
Correction (r 2 ) 
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15a Unfocused SAS Image without Range Correction (r 2 ) 
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1 5b Unfocused SAS Image with Range Correction (r 2 ) 

Figure 15 Unfocused SAS Images of Sphere with and without range 
Correction (r 2 ) 

Figure (15a) shows the strong cross talk between the transducer and the microphone that 
occurs around range bin 5. Range bins 5 to 20 exhibit the effect of the die-out of this cross 
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talk, which can also be seen in Figure (14a). The sphere is located at about range bin 30. The 
effect of the r range correction reduces the amplitude of the cross talk and made the effect 
of the die-out look almost range independent. The ambient noise at further range bins are 
however amplified by the range correction. 

In the case of focused SAS processing, the spherical spreading of the transmit 
pulse is largely recovered by the summation of signals over the 3 dB beamwidth envelope, 
which collects about 87% of the transmitted power. Hence, for SAS processing, the range 
correction only needs to compensate for the divergence from the transmit path. This is shown 
by the images in Figures (16a) and (16b). 
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16a Focused SAS Image without Range Correction (r 1 ) 
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16b Focused SAS Image with Range Correction (r 1 ) 

Figure 16 Focused SAS Images of Sphere with and without Range 
Correction (r 1 ) 
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3. Compressed Signal 

For pulse compression, the received signal was correlated with the transmitted 
pulse. For a given azimuth position, the received signal had 584 data elements while the 
transmitted pulse, also known as the reference pulse, had 65 data elements. Using the 
convolution function in Matlab, the convoluted data had a total of 648 data elements. The 
584 element pulse compressed signal actually started on element 33 and extended to element 
615. This is demonstrated by correlating the reference pulse with itself as shown in Figures 
(17a) and (17b). 
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17b Correlation of Reference Pulse with Itself 
Figure 17 Signal Correlation for Pulse Compression 
The correlation of the reference pulse with itself produces a data array with a total of 129 
data elements. The 65 data element pulse compressed signal was actually from data 
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elements number 33 to 97 as shown in Figure (17b). An error in selecting the starting point 
of the pulse compressed signal will produce an image that is shifted either forward or 
backward in range. For example, in the case of the sphere target and reference pulse of 65 
data elements, choosing the pulse compressed data array to be data element 1 to 584 will 
shift the image by 32 range bins further away. The SAS images of the sphere with 
compressed data array for the different combinations described above are shown in Figure 
(18). 



70 

60 

50 

<D 

[?40 
* 30 
20 

20 40 60 80 100 120 

Azimuth 

18a Chirped-Focused SAS Image of a Sphere starting at Data Elements 33 
to 616. 
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18b Chirped-Focused SAS Image of a Sphere starting at Data Elements 1 
to 584. 
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18c Chirped-Focused SAS Image of a Sphere starting at Data Elements 65 
to 648 

Figure 18 Chirped-Focused SAS Images of Sphere for Different Compressed 
Data Array 
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Besides introducing an incorrect range, the image quality also suffered. This was especially 
obvious from Figure (18c) where the strong return at the far range was from the range 
amplification of the noise growing in comparison to the target signal. 

4. Translation Distance 

For each data run, travel time of the trolley was measured using a stop watch. 
The translation distance that the transducer had travelled was measured using the ruler 
mounted on the track. For 1024 azimuth elements sampled at 18 elements per second, the 
data collection should take about 56.9 seconds. However, the time taken for one data run 
was about 63 seconds. The distance travelled was found to be 9 cm longer than that 
computed for the time to capture 1024 azimuth elements by the A-D converter. Indeed, the 
image improved quite significantly when the computed translation distance was used rather 
than the measured translation distance. Figure (19) shows the chirped-focused SAS image 
of the sphere using the measured distance while Figure (20) show the image using the 
computed distance. 
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Figure 19 Chirped-Focused SAS Image using Measured Distance (0.99 m) 
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Figure 20 Chirped-Focused SAS Image using Computed Distance (0.9 m) 
The reason for the difference was not obvious to the author. What seems to 
be clear is that the computed distance should be used as this is derived from the rate at which 
data is actually captured by the computer. A possible cause could be delays in the start up or 
ending of the data acquisition process. 
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With the above mentioned adjustments, the performance of the different 
algorithms was evaluated. Figure (21) shows the raw, real-aperture image with five cycle, 
sinusoidal pulse transmission at 25 kHz. The presence of a target can be guessed. The 
hyperbolic arc on both sides of the target can be clearly seen. 




Azimuth 

Figure 21 Raw, Real -aperture Image of a Sphere at 25 kHz 
Figures (22) and (23) show the SAS images at 25 kHz for the unfocused and 
focused cases. The unfocused SAS image is a strong confirmation that there was indeed a 
target present. The range resolution was about 5.5 cm and azimuth resolution was about 4.7 
cm. The focused SAS image pinpointed the target with very good resolution. With the 
display scale as shown, the azimuth resolution was indeed 0.8 cm or better while the range 
resolution was 2.7 cm or better. 
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Figure 22 Unfocused SAS Image of a Sphere at 25 kHz 
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Figure 23 Focused SAS Image of a Sphere at 25 kHz 
The range resolution attained with a focused SAS and a five cycle, sinusoidal 
pulse at a fixed frequency of 25 kHz leaves room for improvement. Indeed, the focused 
SAS with up-chirp transmission of frequency range 16 kHz to 32 kHz improved the range 
resolution to 1 .4 cm or better. This is shown in Figure (24). 
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Figure 24 Chirped-Focused SAS Image of a Sphere 
Figure (24) show what appears to be ‘wings’ on the image when pulse 
compression was used in conjunction with the SAS azimuth processing. Apparently, when 
the hyperbolic summation for each pixel in the image included the target, the sum did not 
average to zero, as it should. The peak signal to rms average background in Figure (25) is 35 
dB. 

Besides the improvement in range and azimuth resolutions, there were also 
significant improvements in the signal-to-noise ratios. However, the price of a more complex 
processing algorithm is in the processing time for each image. Table (1) shows the 
approximate signal-to-noise ratios of the sphere image using the various techniques and the 
respective processing times using a 1 00 MHz Intel Pentium PC operating under Windows 
95. 
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PARAMETER 


PROCESSING ALGORITHM 


Real-Aperture 


Unfocused 

SAS 


Focused SAS 


Focused SAS 
with Pulse 
Compression 


Signal-to- 
Noise Ratio 
(dB) 


7 


15 


26 


35 


Processing 
Time (sec) 


9 


35 


121 


156 



Table 1 Comparison of SNR and Processing Times in Seconds 



B. N-TARGET 

The N-target consisted of sixteen 9- volt batteries arranged in a letter N pattern. The 
batteries were all lying on their longer side with the exception of the center battery which 
stood on its end. Each battery was 4.5 cm by 2.5 cm which was approximately three 
wavelengths by two wavelengths. The images for real-aperture, unfocused and focused SAS, 
and focused SAS with pulse compression are shown in Figure (25). 
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25a Real-aperture Image of N-target 



25b Unfocused SAS Image of N-target 
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25c Focused SAS Image of N-target 25d Chirped-Focused SAS Image of N-target 

Figure 25 Images of N-target 

The real-aperture image reveals a clue that there is a target but not what it looks like. The 
unfocused image while showing that it is shaped like an N stops short on the number of 
components making up the target. The focused images with and without pulse compression 
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show that the target is shaped like a letter N and possibly made up of about 1 6 components. 
Each component was shaped like a rectangle. In addition, the center component is the most 
pronounced and seems to be oriented differently from the other components. Indeed, the test 
against this complex target suggests that at least a classification capability has been achieved. 
The subsequent tests evaluated the algorithms against targets more like those found in an 
operational mine countermeasure environment for littoral areas. 

C. ROCKS 

Two porous rocks of about 14 cm by 10 cm were selected as targets since they 
represent typical sea bottom objects. In addition, the rocks provided a quasi diffuse target that 
should show numerous scattering centers needed to create an extended image. The rocks 
were similar to volcanic rocks with holes that were about a wavelength in diameter. The 
images below are for two such rocks positioned closely to each other. Figure (26) show a 
picture of the rocks as positioned in the experiment. A 15 cm ruler is placed in the 
foreground to provide visual appreciation of the size of the rocks. Figure (27) shows the 
images of the two rocks. 



52 




Figure 26 Picture of 2 Porous Rocks with a 15 cm Ruler 
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27a Real-aperture Image of Two Rocks 27b Unfocused SAS Image of Two Rocks 
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27c Focused SAS Image of Two rocks 27d Chirped-Focused SAS Image of Two rocks 

Figure 27 Images of Two Porous Rocks 

Again, the real-aperture image reveals very little detail of the target although it is evident that 
something is there. The unfocused image also gives little clue on the shape of the target. The 
focused images show that there is some structure to the target and it seems that there are two 
targets close to each other. The chirped-focused image shows pretty clearly that there are two 
targets and they have different shape and orientation. 

These targets are analogous to a diffuse surface viewed optically rather than a 
specular mirror surface. What is apparent in the latter two images is some semblance of the 
granular feature due to speckle. This granular characteristic made it difficult to see the shape 
of the focused SAS image. Pulse compression clearly revealed more detail about the size and 
shape of the targets. 
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D. 



MINES 



Perhaps the ultimate test of the processing algorithms was with mines as they are the 
targets that a sonar system would want to detect, classify and identify, and to do so with a 
high level of confidence. Two mines, the TM46 anti-tank mine (ATM) and the PDM-1 anti- 
invasion mine (AIM), were used in the test. These mines are on loan from Naval Coastal 
Systems Station, Naval Surface Warfare Center, Panama City. The pictures of the ATM and 
AIM are shown in Figure (28). 




28a Anti-tank Mine 28b Anti-invasion Mine 

Figure 28 Picture of the Mines 

The images of the two mines are shown in Figures (29) and (30). 
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29a Real-aperture Image of an 
Anti-tank Mine 
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29b Unfocused SAS Image of an 
Anti-tank Mine 
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29c Focused SAS Image of an 
Anti-tank Mine 



29d Chirped-Focused SAS Image of an 
Anti-tank Mine 



Figure 29 Images of a Russian, TM-46 Anti-tank Mine 
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30a Real-aperture Image of an 
Anti-invasion Mine 
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30b Unfocused SAS Image of an 
Anti-invasion Mine 
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30c Focused SAS Image of an 30d Chirped-Focused SAS Image of an 

Anti-invasion Mine Anti-invasion Mine 

Figure 30 Images of a Russian PDM-1 Anti-invasion Mine 
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The images of the mines collected under very high signal to noise ratio conditions are not 
sufficient to identify the targets. Both types of mines are effectively specular to the 
processing algorithms. The target appears as a few bright glints, which do not provide 
sufficient information for classification or identification. The specular characteristic is just 
the opposite of the rough, diffuse, rock targets shown in Figure (27d). 



58 



V. RECOMMENDATION & CONCLUSION 



This thesis investigated a laboratory synthetic aperture sonar designed to test the 
algorithms and techniques needed to detect, classify and identify minelike objects. The FM 
chirped processing algorithm developed in this thesis achieved a range and azimuth 
resolution of about 1 cm. The 2:1 frequency spread of the FM chirp increased the signal to 
noise ratio by 20 dB compared to an unfocused synthetic aperture system. The processed 
images clearly show detection of targets with a high degree of certainty. However, the ability 
to classify and identify mines and rocks is degraded by speckle and glint effects. 

These results reveal at least two important operational implications for sonar 
applications. One is that the threat targets must be analysed and their sonar images collected. 
In this way, a threat library can be used to aid in the classification and identification of 
targets. Another implication is the importance of operator training. While a stealthy mine 
presents scattered ‘point’ information, this information can be utilized to classify the target 
if the operators are trained to read sonar images which are very different from photographs. 
Indeed, with some prior information, both the anti-tank and anti-invasion mines may be 
classified based on the processed images with a high level of confidence. 

Notwithstanding the above, there is considerable need for improvement. The 
processing algorithms have so far been tested only with propagation data collected in air 
medium and in an ideal laboratory environment, rather ideal situations. There is a need to 
evaluate the performance with data from propagation in water. When collecting data in 
water, the lateral movement of the aperture due to underwater wave motion will reduce the 
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lateral coherence length of the acoustic field, which will degrade the image. In this case, 
motion compensation will be necessary as any lateral movement of the aperture will have 
impact on the phase of the signal. The processing has been carried out off-line after the data 
has been collected. The next challenge would be to collect and process the data using the 
same computer in real-time. 

This thesis found that the high resolution algorithms improved the detection of 
targets and image quality with high signal to noise ratio. However, resolution alone was not 
sufficient to classify and identify minelike targets in complex backgrounds. Resolution of 
this problem will require a different approach such as utilizing adaptive acoustic daylight to 
avoid the speckle and glint problems inherent with even a 2:1 FM chirped or ultra-short, 
mono-pulse synthetic aperture sonar. Hence, to achieve a classification and identification 
capability, a completely different approach to acoustic illumination and signal processing is 
needed. 
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APPENDIX A. DATA SAMPLING PROGRAM 



On Testpoint version 1.1, this program samples 1 024 by 584 data elements at a rate of 1 00 



kHz. 



1 ) 

2 ) 

3) 

4) 

5) 

6 ) 

7) 

8 ) 
9) 



Open 

Set 

Set 

Linear series 
Acquire A/D 
Output to 
End 
Close 
Set 



File 1 

A/D trigger digital A/D 1 
Indicator to 1 

Loop 1 from 1 to 1 024, step by 1 

A/D 1 # samples = 584, rate = 100000 Hz, channel(s) = 0 

File 1 with A/Dl, term. - none 

Loop 1 
File 1 

Indicator to 0 
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APPENDIX B. DATA PROCESSING PROGRAMS 



1. UNFOCUSED SAS 

%This is the unfocused SAS algorithm. The approach is to carry out 
%a coherent summing in the azimuth axis before a moving average. 
%27 Feb 1996 
%By LIM CHIN HU AT 

rb=584; 

azb=1024; 

tO=clock; 

fid=fopen('a:spheref.dat'); 

a=fread(fid,'intl6'); 

fclose(fid); 

b=reshape(a,rb,azb); 
clear a; 

r=rand(size(b)); 

b=r*800+b; 

nas=8; 

offset=mean(mean(b(500 : 5 80, 1 0 : 30))); 
b=b-offset; 

%Compressing the received matrix for faster printing 

cazb=azb/nas; 

crb=rb/nas; 

rcv=zeros(rb,cazb) ; 

rcvO=zeros(crb,azb); 

rcvO 1 =zeros(crb,cazb); 

for n=l:crb 
rcv0(n, : )=b(n * nas, :) ; 
end 

for n=l:cazb 
rcv01(:,n)=rcv0(:,n*nas); 
end 
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% Coherent summing over nas azimuth samples 
for n=0:cazb-l 
l=n*nas; 
for m=l:nas 

rcv(: ,n+ 1 )=rcv( : ,n+ 1 )+b(:,l+m); 
end 
end 



% Carrying out a moving average after coherent summing 

rev 1 =zeros(rb,cazb+nas); 
for p=0:nas-l 
for q=l:cazb 

rev 1 ( : ,p+q)=rcv 1 ( : ,p+q)+rcv(: ,q) ; 
end 
end 

% Carrying out magnification of target signal by range (r A 2) component 

rcv2=zeros(rb,cazb+nas); 

for n=l:rb 

rcv2(n,:)=abs(rcvl(n,:))*n*n; 

end 

%Compressing range component for faster printing 
rcv2 1 =zeros(rb/nas,cazb+nas); 
for n=l:rb/nas 
rcv2 1 (n, :)=rcv2(n*nas, :); 
end 

rcv6=rot90(rcv2 1 ); 
rcv7=rot90(rcv0 1 ); 

%Printing out Received Signal 

figure(l) 

subplot(221) 

title('Received Signal') 

ylabel('Range') 

xlabel(' Amplitude') 

plot(rcvOl) 

subplot(222) 

plot(rcv7) 

xlabel(' Azimuth') 
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ylabel(' Amplitude') 

figure(2) 

subplot(221) 

colormap(flipud(gray(256))) 

h=pcolor(rcv01) 

set(h,'EdgeColorVnone') 

set(h,'FaceColor','interp') 

title ('Real-Aperture Image') 

ylabel('Range') 

xlabel(' Azimuth') 

%Printing out Unfocused SAS Image 

figure(3) 

subplot(221) 

title('CH before MA') 

xlabel('Range') 

ylabel(' Amplitude') 

plot(rcv21) 

subplot(222) 

plot(rcv6) 

xlabel(' Azimuth') 

ylabel(' Amplitude') 

figure(4) 

subplot(221) 

colormap(flipud(gray(256))) 

h=pcolor(rcv21) 

set(h,'EdgeColor','none') 

set(h,’FaceColor','interp') 

title ('Unfocused Image') 

ylabel('Range') 

xlabel(' Azimuth') 



etime(clock,tO) 



2 . 



FOCUSED SAS 



%This program is a fully focused SAS algorithm without interpolation between sample 
elements. 

%27 Feb 1996 

%By LIM CHIN HU AT 

rb=584; 

azb=1024; 

fid=fopen('a:spheref.dat'); 

a=fread(fid,’intl 6'); 
fclose(fid); 

tO=clock; 

bc=reshape(a,rb,azb); 

%Correcting for bias in the data 

offset=mean(mean(bc(500:550, 10:30))); %A known subset of the matrix where there is no 
signal 

b=bc-offset; 

% Coherent summing over nas azimuth samples 
nas=8; 

cazb=azb/nas; 

rcv=zeros(rb,cazb); 

for n=0:cazb-l 
l=n*nas; 
for m=l:nas 

rcv( : ,n+ 1 )=rc v(: ,n+ 1 )+b( : ,1+m) ; 
end 
end 

%Creating larger matrix with zero pegging on both sides required for the SAS processing 

c=zeros(rb,cazb+ 1 22); 

d=zeros(rb,cazb); 

c=[zeros(rb,6 1 ),rc v,zeros(rb,6 1 )] ; 

m=l :rb; 
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%mmax=ceil(0.1042*m); %Az correction for 0.85m and 28 deg bw 
%mmax=ceil(0.090*m); %Azimuth correction for 0.984m run 
mmax=ceil(0.08855*m); %Azimuth correction for 1.0m run 

for m=l:rb-l 
rbm=rb-m; 
for p=l :mmax(m) 

% R=7.52*p*p/m; %Range correction for 0.85m run 

% R=10.0778*p*p/m; % Range correction for 0.984m run 

R=10.408*p*p/m; % Range correction for 1.0m run 
if R<=rbm %interpolate for Rp not an integer number 
mml=p; 
end 
end 



q=l:mml; 

% R=7.52.*q.*q./m; %0.85m run 

% R=10.0778.*q.*q./m; %0.984m run 
R=1 0.408. *q.*q./m; %1.0m run 



q=l:mml; 

x2=round(R(q)); 

for n=l:cazb 
nl=61+n; 
q=l:mml; 

b 1 =diag(c(x2(q)+m,n 1 +q)); 
q=l:mml; 

b2=diag(c(x2(q)+m,n 1 -q)); 
d(m,n)=sum(b 1 )+sum(b2); 
end 
end 

% Carrying out signal magification of target signal by range component 
rev 1 =zeros(rb,cazb); 
for n=l:rb 

rcvl(n,:)=abs(d(n,:))*n; 

end 

%Summing over range elements for faster image printout 
crb=floor(rb/ nas) ; 
rcv4=zeros(crb,cazb); 
for n=0:crb-l 
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3. 



CHIRP REFERENCE PULSE 



%This program set up the reference signal that is needed for pulse compression. 

%27 Feb 1996 

%By LIM CHIN HU AT 

rb=584; %Define the number of range bins 
azb=1024; %Define the number of azimuth bins 

fid=fopen('c:\temp\chirp.ref); 
a=ffead(fid,'intl 6'); 
fclose(fid); 

bc=reshape(a,rb,azb); 
clear a; 

%Correcting the data for bias in the measurement 

offset=mean(mean(bc(5 00 : 5 5 0, 1 0 : 3 0))) ; 
b=bc-offset; 

% Averaging over the azimuth axis 

rcv=zeros(rb); 

c=rot90(b); 

rcv=mean(c); 

%Getting the signal for correlation 

refl=rcv(2 19:283); %65 range bins of signal for correlation/convolution 
ref=fliplr(refl); %Needs to flip for correlation 

%Plotting reference pulse and convolution on itself 
figure (1); 
subplot(221); 
plot (ref); 

title('Reference Pulse') 
xlabel(' Azimuth’) 
ylabel('Amplitude') 

subplot(222) 
plot(abs(conv(refl ,ref))) 
title('Convoluted Signal') 
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xlabel('Azimuth') 

ylabel('Amplitude') 
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4. 



CHIRPED-FOCUSED SAS 



%This program is a fully focused SAS algorithm with pulse compression. 

%27 Feb 1996 
%By LIM CHIN HU AT 

rb=584; 

azb=1024; 

%load chirpref.m %Loading the reference file for the correlation pulse named as REF 

fid=fopen('a:spherec.dat'); 

a=fread(fid,'intl 6'); 
fclose(fid); 

tO=clock; 

%Carrying out convolution of the data with the reference signal 
c=conv(a,ref); 

cl=c(33:598048); %Having only the necessary number of elements 

bc=reshape(c 1 ,rb,azb); 

%Correcting for data bias 
offset=mean(mean(bc(500:550, 1 0:30))); 
b=bc-offset; 

% Coherent summing over nas azimuth samples 
nas=8; 

cazb=azb/nas; 

rcv=zeros(rb,cazb); 

for n=0:cazb-l 
l=n*nas; 
for m=l:nas 

rc v( : ,n+ 1 )=rcv( : ,n+ 1 )+b(: ,1+m); 
end 
end 

%Creating larger matrix require for SAS processing 
c=zeros(rb,cazb+l 22); 
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d=zeros(rb,cazb); 

c=[zeros(rb,61),rcv,zeros(rb,61)]; %Quick pegging 
m=l:rb; 

mmax=ceil(0.1042*m); %Az correction for 0.85m and 22 deg beamwidth 
%mmax=ceil(0.093*m); %Approx for 0.95m run 
%mmax=ceil(0.090*m); %Azimuth correction for 0.984m run 
%mm ax=ceil(0. 08855 *m); %Azimuth correction for 1.0m run 

for m=l :rb-l 
rbm=rb-m; 
for p=l:mmax(m) 

R=7.52*p*p/m; %Range correction for 0.85m run 
% R=9.4*p*p/m; %Approx for 0.95m 

% R=1 0.0778 *p*p/m; % Range correction for 0.984m run 

% R=10.408*p*p/m; % Range correction for 1.0m run 

if R<=rbm %interpolate for Rp not an integer number 
mml=p; 
end 
end 



q=l:mml; 

R=7.52.*q.*q./m; %0.85m run 
% R=9.4.*q.*q./m; %Approx for 0.95 

% R=10.0778.*q.*q./m; %0.984m run 
% R=10.408.*q.*q./m; %1.0m run 



q=l:mml; 

x2=round(R(q)); 

for n=l :cazb 
nl=61+n; 
q=l :mml; 

b 1 =diag(c(x2(q)+m,n 1 +q)); 
q=l:mml; 

b2=diag(c(x2(q)+m,n 1 -q)); 
d(m,n)=sum(b 1 )+sum(b2); 
end 
end 

%Carrying out signal magnification of target signal by range component 
rev 1 =zeros(rb,c azb) ; 



72 



for n=l :rb 

rev 1 (n, : )=abs(d(n, : )) * n; 
end 

crb=floor(rb/nas) ; 
rcv2=zeros(crb ,cazb) ; 
for n=0:crb-l 
l=n*nas; 
for m=l:nas 

rcv2(n+ 1 ,:)=rcv2(n+ 1 , :)+rcv 1 (1+m, :) ; 
end 
end 

rcv3 =rot90(rc v2) ; 

figure(l) 

subplot(221) 

plot(rcv2) 

title ('Chirp-Focused Image') 

xlabel(' Azimuth') 

ylabel('Range') 

subplot(222) 

plot(rcv3) 

xlabel(' Azimuth') 

ylabel('Range') 

figure(2) 

subplot(222) 

contour(rcv2) 

title('Chirped-Focused Image') 
xlabel(' Azimuth') 
ylabel('Range') 
subplot(22 1 ) 

colormap(flipud(gray(256))) 

h=pcolor(rcv2) 

title('Focused Image') 

xlabel('Azimuth’) 

ylabel('Range') 

set(h,'EdgeColor','none') 

set(h,'FaceColor','interp') 

etime(clock,tO) 
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APPENDIX C. FM CHIRP PROGRAM 



This FM chirp program is downloaded into volatile memory of the HP33120A waveform 
generator. 



10 IChirp Waveform Generator 

20 OPTION BASE 1 

30 INTEGER Waveform(4000) 

40 INTEGER I,N 

50 ! 

60 ASSIGN @Fgen TO 7 1 0 

70 OUTPUT @Fgen; “*RST” 

80 ! 

90 N=4000 

100 Numcycles=10. ! Number of cycles of start frequency 

110 Scale=Numcycles/N 

120 Fchirp=4./N ! End frequency multiplier 

130 ! 

140 FOR 1=1 TON 

150 F=I 

160 Waveform(I)=2047*SIN(2*PI*Scale*(F+Fchirp*F*F)) 

170 NEXT I 

180 ! 

1 90 ! Download data into volatile memory 

200 ! 

210 DISP “Downloading Arb Waveform” 

220 OUTPUT @Fgen; “DATA:DAC VOLATILE ”;Waveform(*) 

230 DISP “Download complete” 

240 ! 

250 OUTPUT @Fgen;”DATA: COPY CHIRP, VOLATILE” ! Copy ARB to non- volatile 
260 OUTPUT @F gen;”FUNC : U SER CHIRP” 

270 OUTPUT @Fgen;”FUNC:SHAP USER” 

280 OUTPUT @Fgen;”OUTP:LOAD 50" ! Output termination 50 ohms 

290 OUTPUT @Fgen;”FREQ 1 000; VOLT 5" ! Output frequency is 1kHz, 5V 

300 ! 

310 PAUSE 

320 END 
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